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When a low dimensional chaotic attractor is embedded in a three dimensional space its topological 
properties are embedding-dependent. We show that there are just three topological properties that 
depend on the embedding: parity, global torsion, and knot type. We discuss how they can change 
with the embedding. Finally, we show that the mechanism that is responsible for creating chaotic 
behavior is an invariant of all embeddings. These results apply only to chaotic attractors of genus 
one, which covers the majority of cases in which experimental data have been subjected to topological 
analysis. This means that the conclusions drawn from previous analyses, for example that the 
mechanism generating chaotic behavior is a Smale horseshoe mechanism, a reverse horseshoe, a 
gateau roule, an 5-template branched manifold, . . . , are not artifacts of the embedding chosen for 
the analysis. 



I. INTRODUCTION 

Chaotic time series have been generated by a large 
number of experiments. Typically a scalar time series 
is available, and a chaotic attractor must be generated 
from the scalar time series using some embedding proce- 
dure. The algorithm of choice is the time delay embed- 
ding P, [3, although differential embeddings, Hilbert 
transform embeddings, and singular value decomposition 
embeddings have also been used [1, Q • 

The properties of embedded chaotic attractors have 
been analyzed along three distinct mathematical lines: 
geometric, dynamical, and topological. Geometric analy- 
ses involve computing various fractal dimensions fo'] . Dy- 
namical analyses involve computing Lyapunov exponents 
and the average Lyapunov dimension Q. Topological 
analyses concentrate on the global topological properties 
of an attractor by studying how stretching and squeezing 
mechanisms organize the unstable periodic orbits embed- 
ded in the attractor [1 S i, i, . 

In all approaches, it is assumed that the embedding 
adopted creates a diffeomorphism between the under- 
lying (invisible) experimental attractor that generates 
the data and the embedded, or reconstructed, chaotic 
attractor 0, [^. Since the geometric and dynamical 
measures (dimensions and exponents) are invariant un- 
der diffeomorphisms, in principle these real numbers are 
embedding-independent. In practice they are difficult 
to compute, and become increasingly difficult to com- 
pute as the length of the time series and/or the signal 
to noise ratio decreases. Further, there is no indepen- 
dent way to compute errors for the estimates of these 
real numbers. It was even shown in that in some ex- 
perimental data sets estimates of the fractal dimensions 
are diffeomorphism-dependent. Spurious Lyapunov ex- 
ponents occur when the embedding dimension exceeds 
the dimension of the dynamical system. This has been 
addressed in [l^ [isj but remains an open problem. 



By constrast, topological analyses on three dimen- 
sional embeddings have been carried out with relatively 
short experimental data sets and are robust against noise. 
In addition, this analysis method is overdetermined. The 
stretching and squeezing mechanism creating the embed- 
ded chaotic attractor can be determined from a small 
number of unstable periodic orbits and used to pre- 
dict the topological organization of all remaining or- 
bits. These predictions (linking numbers, relative rota- 
tion rates) either agree or do not agree with those for 
orbits extracted from the embedded chaotic attractor. 
In the latter case the model describing stretching and 
squeezing must be rejected. 

What has never been satisfactorily understood is the 
relation between the topological properties of the under- 
lying (invisible) experimental attractor that generates 
the data and the chaotic attractor that has been con- 
structed through an embedding of the data. We illus- 
trate this difficulty with two examples. (1) The Lorenz 
attractor [ij] is described by variables {x{t),y{t), z{t)). 
One three-dimensional embedding is the obvious one: 
{Xi,X2,X^) = {x,y,z). The chaotic attractor in this 
embedding is invariant under rotations: (Xi, X2, X^) — > 
(— Xi, — X2, X3). On the other hand, if a single vari- 
able is observed (either x{t) or y{t)) [l^ the chaotic 
attractor created through the differential embedding 
(Xi, X2, X:i) — (x,x,x) will exhibit inversion symme- 
try {Xi,X2,X3) ^ (-Xi,-X2,-X3) MM- (2) The 
chaotic behavior of Benard-Marangoni fluid convection 
in a square cell [18\ was modeled by a periodically driven 
Takens-Bogdanov nonlinear oscillator [l^. This model 
was studied using a time delay mapping of the form 
{Xi,X2,X3) = {x{t),x{t),x{t - t)) [20]. For a range of 
values of the time delay r, ri < r < T2, the image of the 
data under this mapping exhibits self intersections and 
the mapping is therefore not an embedding (technically, 
it is an immersion) [20, 21]. For r < ti and T2 < t the 
mapping is an embedding. The topological organization 
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of the unstable periodic orbits under the two embeddings 
is different, so the global topological structure of the two 
embedded attractors is not equivalent [2l| . 

This discussion brings us to the crucial question: 
When topological information about a chaotic attrac- 
tor is determined from a three-dimensional embedding 
of the chaotic attractor, what part of that information 
is embedding-dependent and what part is embedding- 
independent? In this work we answer this question for 
a large class of chaotic attractors. These consist of all 
chaotic attractors of "genus-one" type [13, : their nat- 
ural phase space is equivalent to a torus. This includes 
the chaotic attractor discussed in the example (2) above. 
The answer is that the "mechanism" (defined in Sec. |V] 
below) is independent of embedding. Further, the topo- 
logical organization of all periodic orbits in the attractor 
can differ in a very limited number of ways (parity, global 
torsion, and knot type, see Sec. IIVI below). This crucial 
question remains open for chaotic attractors whose nat- 
ural phase space is a torus of genus g {g > 1). This in- 
cludes the Lorenz attractor as well as many other chaotic 
attractors [13, [13] ■ It also remains open for all higher- 
dimensional (hyper-)chaotic attractors. 

II. ASSUMPTIONS 

We make the following assumptions: 

1. A deterministic process (e.g., laser equations, 
Navier-Stokes equations) acts to generate an exper- 
imental chaotic attractor that is three-dimensional. 
A single variable (e.g., laser intensity, fluid surface 
height) is measured. 

2. At least one embedding of this scalar time series 
in can be constructed. This embedding creates 
a diffeomorphism between the original experimen- 
tal chaotic attractor and the embedded or "recon- 
structed" chaotic attractor. 

3. The embedded chaotic attractor is of genus-one 
type: that is, it can be enclosed in a genus-one 
bounding torus [H, [1^ . 

Some remarks about these assumptions are in order. 
We assume in (1) that there is an experimental chaotic 
attractor and that it is three dimensional. By three- 
dimensional we mean explicitly that there is a three di- 
mensional manifold in the phase space that contains the 
attractor. We require this assumption on dimension be- 
cause, at the present time, topological analysis methods 
based on templates are only applicable to three dimen- 
sional chaotic attractors, that is, those that that exist in 
three-dimensional manifolds. The assumption that the 
deterministic process generates a low-dimensional attrac- 
tor is also strong: the Navier-Stokes and the full laser 
equations are partial differential equations rather than 
sets of ordinary differential equations, and act in Hilbert 
spaces rather than finite dimensional phase spaces j25| . 



Assumption (2) is necessary because the Whitney em- 
bedding theorem [lij and its dynamical variants [2, 3] only 
guarantee that the three-dimensional manifold contain- 
ing the chaotic attractor can be embedded into a space of 
sufficiently high dimension (6 = 2x3), but do not ensure 
that it can be done into a three-dimensional phase space. 
In practice, whether this assumption holds can be tested 
a posteriori by verifying that the topological invariants 
measured are consistent with a single two-dimensional 
branched manifold. The diffeomorphism property that is 
assumed of the mapping is the standard assumption for 
all approaches to analysis of embedded data [3,1^. 

Assumption (3) is crucial for our result. It allows us to 
reduce the problem of the inequivalence of embeddings 
of chaotic attractors to the problem, already solved [2q . 
of the equivalence classes of diffeomorphisms of the solid 
torus into the three-dimensional euclidean space . In 
the higher genus case (e.g., Lorenz attractor) the spec- 
trum of inequivalent diffeomorphisms (embeddings) of 
these attractors is related to the spectrum of inequivalent 
diffeomorphisms of the higher genus tori to themselves, 
which remains to be studied. 



III. PRELIMINARY REMARKS 

We begin by recalling that diffeomorphisms map pe- 
riodic orbits to periodic orbits. If x(t) is a point on a 
periodic orbit so that x(t + T) = x(i), then under a dif- 
feomorphism that takes x — > y, y(t) = y{t + T). This 
means that periodic orbits are neither created nor anni- 
hilated by diffeomorphisms. In particular, the spectrum 
of periodic orbits associated with ("in") a chaotic at- 
tractor is an invariant of diffeomorphisms. On the other 
hand their topological organization, as encoded by their 
topological invariants (linking numbers, relative rotation 
rates) could change under diffeomorphism. 

We will describe exactly how these topological invari- 
ants can change under diffeomorphism when the phase 
space containing the chaotic attractor is a torus x S*^ , 
where is a disk in the plane [D^ C R^) and S*^ is 
parameterized by cj), usefully considered as a phase angle 
mod 27r. In this phase space trajectories can be expressed 
in the form {x{t),y{t),(j){t)). This class includes nonau- 
tonomous dynamical systems such as the periodically 
driven Duffing, van der Pol, and Takens-Bogdanov non- 
linear oscillators where and t are linearly related, and 
autonomous dynamical systems whose phase space pro- 
jection (x, x) exhibits a 'hole in the middle' (e.g., Rossler 
system ^ at (a,5,c) = (0.398,2.0,4.0)). It includes 
other autonomous dynamical systems with a hole in the 
middle that is present but obscured by simple projections 
(e.g., Rossler system at (a, 6, c) = (0.398,2.0,13.3)). For 
this class of systems the phase (j) = (j){t) is a monotonic 
function of the time t. This discussion explicitly excludes 
attractors of genus g > 2 with two or more 'holes in the 
middle', such as the Lorenz attractor. 

In the work to follow we seek a discrete enumeration 
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of embeddings, or difTeomorphsims, of strange attractors. 
To achieve this end it is necessary to "mod out" continu- 
ous degrees of freedom associated with diffeomorphisms. 
To do this we introduce the idea of isotopy. Two embed- 
dings /o and /i are isotopic is there is a one parameter 
family of mappings, /(s), with /(O) = /o, /(I) = /i 
and /(s) is an embedding for all s, < s < 1. Such 
a family of embeddings merely deforms the phase space 
smoothly. The topological organization of periodic orbits 
is unchanged under isotopy. For if two orbits intersected 
during the deformation from s = to s = 1 the unique- 
ness theorem would be violated and the mapping /(s) 
(for some s) would not be a diffeomorphism. For this 
reason isotopic mappings are in some sense equivalent. 
The sense is that all topological indices for orbits in a 
strange attractor are the same for all embeddings in the 
same isotopy class. 

Our problem therefore reduces to (1) classifying the set 
of isotopy classes of diffeomorphisms x 5^ — > x 5'"'^, 
(2) classifying the set of isotopy classes of diffeomor- 
phisms X — > i?"^, and (3) determining how topo- 
logical invariants change from one class to another. The 
first two parts of this program are resolved in Sees. IVa 
and IVb. The third part is discussed in Sees. V and VI. 
A more detailed exposition of these points is presented 
in Appendix A. 



IV. EMBEDDINGS OF A TORUS 

Diffeomorphisms of the torus fall into two broad 
classes: intrinsic and extrinsic [1^ [2^. Intrinsic dif- 
feomorphisms are mappings of the torus to itself "as 
seen from the inside." Specifically, they are mappings 



X 



D X S . Extrinsic diffeomorphisms de- 



scribe how the torus sits in R . They are mappings 
X R^. Intrinsic diffeomorphisms are respon- 

sible for two of the three degrees of freedom mentioned 
in the abstract and introduction: parity and global tor- 
sion. Extrinsic diffeomorphisms are responsible for the 
first two and in addition the third: knot type. 



A. Intrinsic Diffeomorphisms. 

These also fall into two classes: those that are isotopic 
to the identity and those that are not. 

Isotopic to the Identity. Diffeomorphisms that are 
isotopic to the identity smoothly deform the phase space. 
Therefore they do not change the topological organiza- 
tion of the periodic orbits in it. Under these diffeomor- 
phisms the topological invariants of periodic orbits re- 
main unchanged. 

Not Isotopic to the Identity. Mappings of the 
torus to itself that are not isotopic to the identity have 
been classified 26]. The idea is as follows. On the two- 
dimensional surface — d{D^ x S^) that is the bound- 
ary of the solid torus it is possible to construct two circles 



that cannot be deformed to a point, as shown in Fig. [T] 
We orient both. The longitude is oriented along the direc- 
tion of the dynamical system flow. The meridian bounds 
a disk that can be used as a Poincare surface of section. 
It is oriented according to the right hand rule. Up to 
isotopy (the class of diffeomorphisms considered in the 
preceeding paragraph) the inequivalent diffeomorphisms 
of the torus to itself are classified by their action on the 
longitude and meridian by the matrix j26|] 



1 n 
±1 



(1) 



The integer n describes the number of rotations of the 
longitude about the core (center line) of the torus as the 
phase angle (f) increases from to 27r. The integer ±1 in- 
dicates whether the diffeomorphism preserves or reverses 
the orientation of the meridian. We identify ±1 with 
parity and n with global torsion in Sec. [V] 






(b) 



(c) 



FIG. 1: (a) Two nonisotopic circles are drawn on the surface 
of the solid torus containing a chaotic attractor. The longi- 
tude is oriented along the direction of the flow. The meridian 
is oriented by the right hand rule. The solid torus is mapped 
diffeomorphically to a torus with (b) n = 2 or (c) negative 
parity. 

Remark. The matrices presented in Eq.([T]) are group 
operations. Diffeomorphisms of the torus to itself form a 
group. The subset that is isotopic to the identity forms 
a subgroup that is invariant in the larger group. The 
quotient of these two groups therefore forms a group. 
This group is discrete. It is generated by two operations, 
represented by the matrices 



1 1 
1 



and 



1 

-1 



(2) 



The first describes the generator that produces a uni- 
form rotation along the axis of the torus: ((a; -t- iy), (p) — > 
{{x + iy)e^'^,(j>). The second generator produces the ef- 
fect of looking into a mirror: {x, y, (p) — > {x, —y, 0). This 
coset decomposition says simply that every intrinsic dif- 
feomorphism can be constructed by composing a dif- 
feomorphism isotopic to the identity with one from the 
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discrete group whose matrix representation is given in 
Eq.©. 



B. Extrinsic DifFeomorphisms. 

The mapping of x into B? shown in Fig. 
is called the 'natural embedding' One nat- 

ural embedding of a chaotic attractor with coordinates 
(xi(0), 0:2(0), (^) in £|2 X into is {X {t) ,Y {t) , Z [t)) , 
with t = (j) and X = [R — xi) cos (j)^ Y ^ [R — xi) sin (/>, 
and Z — X2. This is an embedding provided the circle 
is "bigger" than the attractor. Specifically, if the ra- 
dius of the disk containing the attractor is a, so that 
^\ (4') + ^2 < fo'" ^-ll then R > a guarantees that 
no self-intersections occur in the natural embedding. 




FIG. 2: (a) The torus x of Fig. [Ha) is embedded in 
a natural way in R^. In this embedding the core of the torus 
is a circle of radius R. the torus can also be mapped into R* 
with a nonzero framing index /. The framing index is -f 1 in 
the embeddings (b) and (c). 

The circle is the simplest knot in R"^. Other knots 
in i?"^ can be used as central curves for other extrin- 
sic embeddings. The knot K. has coordinates K(0) = 
{Ki{4>),K2{(t)),K3{(l))) withK((?!)) =K{(j> + 2TT). As with 
any smooth space curve [29] this knot has a moving coor- 
dinate system (repere mobile) with orthogonal unit vec- 
tors t{(j)),n{(j)),h{(j)). The section of a chaotic attractor 
in X S-^ at phase angle (j) is lifted into the plane in 
R'^ perpendicular to the tangent vector t{(j>) at K(0) by 



the mapping (xi ((/)), a;2(0), K(0) + xi{(j))n{(t)) + 

a;2(0)b((/)). This mapping is an embedding provided 
there are no self intersections. This is guaranteed pro- 
vided two conditions are satisfied [30l |: 

Local condition: The radius of curvature of K. is every- 
where greater than a. 

Global condition: The curve IC is "big enough." This 
means specifically that all nonzero local minima of 
|K(0i) — K(02)| are larger than 2a. 

An important integer is associated with each knot IC. 
This is its framing index, / [2^. It describes how many 
times the vectors n and b wind around t as the knot 
is traversed. Specifically, it is the gauss linking num- 
ber of two closed curves in R'^. One closed curve is the 
knot itself. The other is obtained by displacing it a small 
distance along the normal vector. Its coordinates are 
given by setting (xi ((/)), 0:2 (</'),</') = (1,0,(/)) in the map- 
ping above. We use this integer in Sec. IVII to describe 
the problems of the delay embeddings of the fluid data 
presented in Sec. HI Example (2) (embedding of Benard- 
Marangoni fluid data) . Embeddings of the torus into R^ 
with framing index / = +1 are shown in Fig. [2jb) and 
(c). 

Remark. As Fig. [2] shows, choice of a knot in R^ 
for the center curve of the embedded torus is indepen- 
dent of the choice of the framing index of the embedded 
torus. The knot type of the center curve is one degree 
of freedom of embeddings of a genus-one torus into R^. 
Two other degrees of freedom, the framing index (which 
is equivalent to global torsion) and parity have already 
been encountered in diffeomorphisms xS^ ^ x S^. 

Remark. It is pedantically more accurate to describe 
extrinsic embeddings as diffeomorphisms xS^ x 
C R^. In the remainder we forgo this mathematical 
precision. 

V. MECHANISMS 

Chaotic attractors in three dimensional spaces are 
characterized by the spectrum and topological organiza- 
tion of their unstable periodic orbits (UPOs) 0, [a]- The 
topological organization of the periodic orbits is summa- 
rized by a knot holder (also called a branched manifold 
or a template) [3l|, |32|- The spectrum of UPOs in the at- 
tractor is a subset of all the orbits contained in the knot 
holder. This subset is specified by a basis set of orbits 
[3^ . The knot holder that describes an embedded chaotic 
attractor is identified by extracting a rather small set of 
orbits from the attractor and determining their topologi- 
cal organization flO]. As a result, knot holders are invari- 
ant under diffeomorphisms isotopic to the identity, since 
they are derived from the topological indices of periodic 
orbits, which do not change under isotopy. Knot holders 
can differ only by the indices that describe the distinct 
equivalence classes of diffeomorphisms. These are: the 
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parity index ±1, the global torsion n, and the knot type 
of the embedding into R^, including the framing index 
/. Further, as the spectrum of UPOs in a chaotic attrac- 
tor is a diffeomorphism invariant, every embedding of a 
chaotic attractor has the same basis set of orbits. 

A knot-holder has as many branches as the number of 
symbols required to uniquely identify the unstable peri- 
odic orbits in the attractor. This number, as well as the 
symbolic name of each periodic orbit, can be determined 
by constructing a generating partition [s^ [ssIIsgIIstIIsII . 
Techniques have also been developed to construct the 
knot-holder without prior knowledge of a symbolic encod- 
ing, by searching directly for the simplest template with 
a set of orbits isotopic to the experimental one H, HI, \^ . 
A generatin g pa rtition can then be constructed from this 
information [si,!!^,!!!. 

A knot-holder has one or more branch lines. Two or 
more branches leave from each branch line ("stretching 
process" ) , and two or more branches meet at each branch 
line ( "squeezing process" ) . Since knot holders are surro- 
gates for chaotic attractors [Sll, [13] , we regard informa- 
tion about which branches leave each branch line and 
which meet at each branch line as describing the mecha- 
nism generating chaos. 

Chaotic attractors in a torus (genus one) possess a sin- 
gle branch line [22], which may be an interval (Rossler 
and Dufhng attractors) or a circle (van der Pol attrac- 
tor). For attractors in x by "mechanism" we mean 
explicitly the order in which branches leave the branch 
line (left to right) or circle (clockwise or counterclock- 
wise) and the order in which the branches are squeezed 
together when they return to the branch line (front to 
back) or circle (inside to outside) In the genus one 
case, mechanism describes how the branch curve (line, 
circle) is folded back into itself in one forward iteration. 
The return flow, from the output side of the branch line 
(lines c to d in Fig. [3lja)) to the input side (lines a to 6 
in Fig. [3Ka)) is assumed to preserve order. The "mech- 
anism" is shown within the dashed rectangle of Fig. [31 
The part of the branched manifold describing the flow 
from 6 to c is the part of the branched manifold that de- 
scribes stretching (the divergence of branches A and B) 
and squeezing (the joining of branches A and B). This 
is the part of the branched manifold describing "mech- 
anism." This knot-holder has only one branch line. We 
have shown four in Fig. [3] to emphasize the various roles 
played by that branch line. 

Knot-holders for chaotic attractors in a genus-one 
torus are classified by a pair of matrices [l,[l,[a[I3]- If n 
symbols are required to label periodic orbits, one matrix 
(the "template matrix" ) is an ri x n matrix and the other 
("array matrix" or "joining matrix") is a 1 x n matrix. 
These two matrices are shown in Fig. 2] for two particu- 
lar knot-holders. One (Fig. [DJb)) is the outside to inside 
scroll template with three branches, which has been ob- 
served i n ( embeddings of) experimental data from lasers 
[Ulillil] and from neurons [soj. The other (Fig. H^a)) 
is the inside-to-outside gateau roule. The diagonal ma- 
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FIG. 3: (a) Knot holder for the Rossler attractor, shown inside 
a torus x 5*^. The flow enters at a, is split at 6, joined at 
c, and "leaves" at d. Periodic boundary conditions identify 
a and d. We also identify b with a and c with d. (b) Knot 
holder for the Lorenz attractor, shown inside a genus-three 
torus. Branches A and B split at a while C and D split at c. 
Branches C and B join at 6 while D and A join at d. In both 
cases the mechanism is shown within the dashed box (a) or 
boxes (b). 



trix elements Ta of the template matrix describe the local 
torsion (measured in units of tt) for branch i. The off- 
diagonal matrix elements T!y = 2 x Link{i,j) are twice 
the linking numbers of the period-one orbits in branches 
i and j. The array matrix describes the order in which 
the branches are glued together at the branch line: the 
smaller the integer entry, the further from the viewer in 
the projection. 

Mechanisms that differ by being mirror images or by 
having integer global torsion are represented by closely 
related matrices. In the opposite parity case, the mirror 
image knot holder has all integer entries with opposite 
signs. In the case of global torsion n, the even integer 
2n is added to all entries in the template matrix. The 
matrices that describe these two variations of the gateau- 
roule mechanism [c.f. Fig. Sfb)] are 



u 



1 

(a) 



1 

(b) 



Branch Matrices Branch Matrices 



2 
12 
2 2 2 

12 0' 





12 

2 2 

2 1 



FIG. 4: Three-branch knot-holder for an (a) inside-to-outside 
and (b) outside-to-inside "jelly roll" mechanism, the template 
matrices and arrays that describe these branched manifolds 
algebraically are also shown. 



Branch Matrices 
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FIG. 5: Three-branch knot-holder for an S-mechanism, along 
with its template matrix and joining array. 
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[0 


-2 


-1] 



Global Torsion n 

2n 2n 2n 
2n 2n + l 2n + 2 (3) 
2n 2n + 2 2n + 2 

[0 2 1] 



Embeddings with nontrivial knot type do little to alter 
the matrices that describe the mechanism that generates 
chaotic behavior [2l[. Nontrivial knot type may change 
parity and add global torsion, depending on the framing 
[iil of the knot (see Sec.|IV|. 

To be explicit, a mechanism that generates chaos re- 
quiring three symbols can be of two types: a scroll mech- 
anism (Fig. [4]) or an "S" mechanism. The template for 
the latter is shown in Fig. \5\ along with its description 
in terms of matrices. 

If one embedding of data reveals a scroll template, all 
embeddings will reveal a scroll mechanism. If on the 
other hand one embedding reveals an S mechanism, ev- 
ery other embedding of these data will also reveal an S 
mechanism. This is true because no transformation in- 
volving sign changes or addition of global torsion [c.f., 
Eq. can change the description given in Fig. [4] to 
the description given in Fig. [S] The mechanism (scroll 



or S) is an invariant of embeddings. 

Similarly, a horseshoe mechanism 7i(n, e) will be de- 
scribed in all embeddings by template matrices 



Horseshoe with 

Parity and 
Global Torsion 



Branch 


1 



Matrices 

2n 2n 
2n 2n + e 

[0 e] 



(4) 



with e = ±1 and n indicating parity and global torsion 
or framing index, respectively. Matrices (jl]) describe all 
possible templates with two branches folding over each 
other. A mechanism identified as a horseshoe in one em- 
bedding is a horseshoe in any embedding. 



VI. TOPOLOGICAL INDICES 



Relative rotation rates are the natural topological in- 



dex for periodic orbits in the torus D 



Linking 



numbers are the natural topological index for periodic 
orbits in R^. 

Assume A and B are two periodic orbits in some em- 
bedding in the torus x S^, and that their relative ro- 
tation rates are Rij{A, B). These fractions are invariant 
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under diffcomorphisms isotopic to the identity. Under 
difFeoniorphisms x ^ x 5^ with global torsion 
n, or with parity —1, that map A ^ A' and B B' 



Global torsion 
Parity = 



n R,j (A', B') = R^j {A, B)+n 



-1 R^,{A',B') = -R,,{A,B) 

Diffeomorphisms x ^ R^ map A A" and 
B B" . For these closed curves in R^ it is possible to 
compute both relative rotation rates and linking num- 
bers. Under the natural embedding :8] (c.f., Fig. ^ 

R,j{A",B") ^ R^J{A,B) 

Under an embedding into with framing index / 

R,j{A",B") ^ R,,{A,B) + f 

In all cases the linking numbers of A" and B" in R^ are 
the sum of their relative rotation rates [1] : 

PA PB 

L(J^\B")^Y.T.R^M"^B") 

i=i j=i 

where pA = Pa" is the period of orbits A and A" , and 
similarly for B. The dependence of the linking number 
of A" and B" on the framing index / is 

Lf{A",B")^LoiA",B") + fpAPB 

The framing index / for embeddings x S*^ — > x 
C R^ can be considered, for all practical purposes, as 
equivalent to the global torsion n for embeddings x 
S^-.D^xS\ 



VII. PERESTROIKAS 

Up to this point the discussion has concentrated on 
embeddings of a single attractor. Usually experiments 
that generate chaotic attractors are carried out over a 
range of control parameter values in an effort to create 
the equivalent of a bifurcation diagram. In this section we 
discuss fixed embeddings of a family of attractors and the 
dual process: families of embeddings of a single attractor. 

The first topological analysis of a family of chaotic 
attractors was carried out in [i^. A single embedding 
was used to analyze many data sets from lasers with sat- 
urable absorbers operated with three different absorbers 
and under various operating conditions. This analysis 
revealed that through all these changes the underlying 
branched manifold never changed: it was only the basis 

Results for an NM R 



set of orbits that changed H, 
laser [46] and a nonlinear vibrating string [47| were the 
same. Subsequently, studies of the periodically driven 
Duf fing oscillator [4^, CO2 lasers with modulated losses 
[3i,|ii, an iVd-doped YAG laser[i3[^n TVd-doped fiber 
laser (43l |4^ , and sensory neurons 50] showed that the 



underlying branched manifold was a "gateau roule" or 
"jelly roll" branched manifold [1, [E], and that under vari- 
ation of the modulation frequency the flow was directed 
to branches of this branched manifold with systematically 
increasing torsion. 

In light of the results presented in the preceeding 
sections, these conclusions are embedding-independent: 
they would have been reached using any embedding. 
First, the variation of torsion with control parameters 
was observed using a fixed embedding, hence is due to 
physical effects and not to the choice of embedding. 
Within a fixed-embedding study, the standard horseshoe 
H{0, 1) is topologically distinct from a "reverse" horse- 
shoe 1) (as observed in Second, the spiral 
structure that globally describes attractors observed at 
different control parameters values would not have been 
affected if an embedding with different knot type, torsion 
and parity had been chosen. 

It is often the case that families of mappings are stud- 
ied in an effort to identify an optimum embedding. The 
method of minimum mutual information [5l[ was devel- 
oped for precisely this reason. The first systematic study 
of the way the topological properties of an embedded at- 
tractor can depend on the embedding, or change with the 
embedding parameters, was carried out in [21|. This is 
the example (2) summarized in the Introduction. Map- 
pings with a delay t < ti provided embeddings, as did 
mappings with t > T2. In both cases, changing the delay 
T by a little had no effect on the topological indices of 
the periodic orbits. In both cases the torus embedded 
in R^ wound around a vertical axis three times before 
closing. In the transition from one regime of embed- 
dings to the other all relative rotation rates changed by 
±2 (depending on whether the time delay t increases or 
decreases). In the interval ti < t < T2 the mapping 
exhibited self intersections, of the type indicated by the 
arrows in Fig. [UJc). The knot type of the embedding into 
jj>3 j-eniained unchanged but its framing in R^ changed. 
Further, the change was by an even integer. This is a 
signature for framing changes caused by change in hand- 
edness of writhe Q . 



VIII. SUMMARY AND CONCLUSIONS 

When a low dimensional chaotic attractor is embed- 
ded in a three dimensional space, its topological prop- 
erties depend on the embedding. We show that, for a 
large class of low dimensional attractors there are three 
topological properties that are embedding-dependent and 
one that is embedding-independent. The embedding- 
dependent properties are: parity, global torsion, and knot 
type. In the latter case (of mappings x 5^ — > R^), 
the framing index is the global torsion. The embedding- 
independent property is the mechanism that acts in phase 
space to create the chaotic attractor. Mechanism is de- 
fined in Fig. [3] in terms of branched manifolds. The 
class of chaotic attractors for which these results hold in- 
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eludes all genus-one attractors: those whose phase space 
is equivalent (diffeomorphic) to a torus x . This 
class includes the Rossler attractor, periodically driven 
two-dimensional nonlinear oscillators such as the Duff- 
ing, van der Pol, and Takens-Bogdanov attractors, and 
most of the experimentally generated chaotic attractors 
that have been studied by topological methods. The prin- 
cipal result is that any single embedding of a three di- 
mensional attractor in this class suffices to determine the 
mechanism that has generated the chaotic data. This 
class does not include the Lorenz attractor and other at- 
tractors with more than one "hole in the middle." 



APPENDIX A: CLASSIFICATION OF 
EMBEDDINGS OF x INTO . 

1. Introduction 

In this appendix, we provide the interested reader with 
more details about how embeddings of genus-one attrac- 
tors can be classified in terms of knot type, torsion and 
parity. 

Assume that two embeddings and ^'2 of a chaotic 
attractor are possible. The simplest case is when ^1 and 
^2 are isotopic: one embedding can be deformed con- 
tinuously into the other. Equivalence of the topological 
properties of the two embeddings then trivially follows 
from the invariance of the topological indices of closed 
curves with respect to smooth deformations that do not 
induce self-intersections. 

When is not isotopic to 5*2, we exploit the assump- 
tion that the original strange attractor can be enclosed in 
a genus-one torus. We first note that a diffeomorphism 
(or homeomorphism) mapping the original attractor to 
a reconstructed attractor is defined on neighborhoods of 
these two strange sets, and can easily be extended to a 
diffeomorphism (or homeomorphism) between solid tori 
contained in these neighborhoods and enclosing the at- 
tractors. 

Since isotopic embeddings are equivalent, determining 
how topological properties of two genus-one embeddings 
of an attractor can differ thus simply amounts to studying 
isotopy classes of embeddings of Z?^ x into K.'^. There 
are two levels in the classification of these isotopy classes, 
because there are two ways in which two embedded solid 
tori can be non-isotopic. The first level is extrinsic and 
deals with how the core of the solid torus is embedded 
in M.^ . When shrunk to their respective cores, two solid 
tori are isotopic if they have the same knot type. The 
second level is intrinsic and deals with how torus bound- 
aries d{D^ X S^) = are mapped to torus boundaries. 
Two embeddings such that the cores of the embedded 
tori have the same knot type can still be non-isotopic if 
the homeomorphism mapping the boundary of one torus 
to the boundary of the other is not isotopic to identity. 

From this classification, we finally conclude that there 
are three degrees of freedom in which two embeddings 



of a genus-one attractor into can differ: knot type 
(extrinsic level), torsion and parity (at both extrinsic and 
intrinsic levels). 



2. Extrinsic level: knot type 

A necessary condition for two embeddings of a mani- 
fold M to be isotopic is that their restrictions to a given 
submanifold M' C M are isotopic [2^. In particular, 
consider the core of the solid torus x S"^, i.e., the sub- 
manifold C — {B} X with a base point B £ D^. Two 
isotopic embeddings of x into K"^ must embed C 
into M'^ isotopically. Since C is an embedding of 5^ into 
X , embeddings of C in R'^ can be classified as em- 
beddings of in M^, i.e., as ordinary knots. Two torus 
cores are thus isotopic if and only if they have the same 
knot type. Conversely, two embedded solid tori whose 
cores are knotted in different ways cannot be isotopic. 




FIG. 6: Two embeddings of x as solid tori in R'^ cannot 
be isotopic if their cores are not isotopic. 

The actual situation (embedding) of tori in R^, as de- 
scribed by the knot type of the torus core, is called the 
extrinsic structure [28|. Every (tame) knot can be used 
as a centerline for a torus that is embedded in M.^ . 

Assume that two embedded solid tori have isotopic 
cores. This allows us to superimpose the boundaries of 
the two solid tori by isotopy deformations. This does not 
imply that the two embeddings are isotopic. However, 
this indicates that we can now study the classification 
of embeddings at an intrinsic level, by considering map- 
pings of the solid torus into itself and forgetting about 
position in M.^ (looking now at the torus from the inside 
rather than from the outside). Thus, knot type of the 
torus core captures all information about isotopy classes 
at the extrinsic level. 
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FIG. 7: The boundaries of two embedded solid tori with iso- 
topic cores can be superimposed by isotopy deformations. 



3. Intrinsic level: global torsion and parity 

We consider now two embeddings 5*1 and of x 
such that the cores of the embedded soHd tori are iso- 
topic and their boundaries are superimposed (Fig. [7|) . A 
necessary condition for the two embeddings to be iso- 
topic is that their restriction to the boundary of the tori 
are isotopic, as with any submanifold. This condition is 
equivalent to requiring that the restriction of (5'i)^^5'2 
to the boundary of the first torus is isotopic to identity. 
It is also a sufficient condition because a homeomorphism 
of the torus that has its restriction isotopic to identity is 
isotopic to identity [2^. Thus, we are left with study- 
ing isotopy classes of homeomorphisms of the boundary 
d{D^ X S^) = T'^ into itself. 

The group of homeomorphisms of this surface to it- 
self modulo isotopically equivalent embeddings is called 
the mapping class group and is equivalent to the mod- 
ular group of 2 X 2 matrices GL(2; Z) = ^ ^ i with 

a, b, c, d integer and ad — be = ±1 [26(|. This group de- 
scribes how closed curves C are mapped to closed 
curves in under the homeomorphism. The description 
is given in terms of the basis set of (two) loops for the 
homotopy group of . These two cycles are the merid- 
ian and the longitude. The meridian can be regarded as 
the small loop that goes around a tire "the short way" 
and the longitude as a long loop that goes around the tire 
"the other way" (c.f., Fig. [ija)). At the topological level 
(homeomorphism) they are more or less equivalent. At 



the level of dynamical systems they are not. The merid- 
ian bounds a disk that lies inside x and can be 
taken as everywhere transverse to the flow that gener- 
ates the strange attractor or its embedding. This disk 
can be chosen as a global Poincare surface of section. 
The longitude can be chosen in the direction of the flow. 
By restricting to the topology underlying the dynamics, 
we investigate the class of inequivalent homeomorphisms 
of the torus boundary into itself. These are described by 
modular group operations of the form 



M = 



1 n 
e 



(Al) 



with e — ±1. We interpret the integer n as the number 
of times the longitude links the core of the solid torus 
X . Dynamically, n is the global torsion of the 
embedding. The sign e = ±1 identifies the parity of the 
torus homeomorphism. 

A point has to be made regarding parity. The mir- 
ror image of an embedding is also an embedding, which 
differs from the original embedding only by orientation. 
In the mirror image of an embedding, all the topological 
invariants are multiplied by —1. Since an embedding and 
its mirror image cannot be isotopic because orientation is 
preserved under isotopy (orientation cannot change with- 
out inducing self-intersections at some stage of the defor- 
mation) , parity has to be taken account when classifying 
embeddings of solid tori into R'^. 

Mirror image transformations act both at the extrinsic 
and intrinsic levels. Their action at the extrinsic level is 
easily incoporated in the the knot type, which is changed 
into its mirror image. At the intrinsic level, parity is 
taken into account through the sign of the lower diagonal 
entry e of the modular transformation (jAip . which is also 
its determinant. 



4. Summary 

Embeddings of genus-one attractors can be classified 
by studying isotopy classes of the cores of the embed- 
ded tori and of their boundaries. There are three degrees 
of freedom by which embeddings of genus-one attractors 
into M'^ can differ: knot type of extrinsic embedding, 
global torsion, parity (handedness). 
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